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ABSTRACT 

Several lines of evidence suggest that fine silicate crystals observed in primi- 
tive meteorite and interplanetary dust particles (IDPs) nucleated in a supersatu- 
rated silicate vapor followed by crystalline growth. We investigated evaporation 
of //m-sized silicate particles heated by a bow shock produced by a planetesi- 
mal orbiting in the gas in the early solar nebula and condensation of crystalline 
silicate from the vapor thus produced. Our numerical simulation of shock-wave 
heating showed that these /xm-sized particles evaporated almost completely when 
the bow shock is strong enough to cause melting of chondrule precursor dust par- 
ticles. We found that the silicate vapor cools very rapidly with expansion into 
the ambient unshocked nebular region; the cooling rate is estimated, for instance, 
to be as high as 2000 K s _1 for a vapor heated by a bow shock associated with 
a planetesimal of radius 1 km. The rapid cooling of the vapor leads to nonequi- 
librium gas-phase condensation of dust at temperatures much lower than those 
expected from the equilibrium condensation. It was found that the condensa- 
tion temperatures are lower by a few hundred K or more than the equilibrium 
temperatures. This explains the results of the recent experimental studies of 
condensation from a silicate vapor that condensation in such large supercooling 
reproduces morphologies similar to those of silicate crystals found in meteorites. 
Our results suggest strongly that the planetesimal bow shock is one of the plau- 
sible sites for formation of not only chondrules but also other cosmic crystals in 
the early solar system. 
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Introduction 



Vapor-solid (VS) growth is a major process for dust formation in the inner region of 
the early solar nebula, where the gas pressure was too low for a liquid phase to exist stably. 
Actually, there are several lines of evidence in interplanetary dust particles (IDPs) and 
primitive meteorites that the dust condensed directly from the vapor phas e as suggested 



for ex ample by enstatite whiskers elongated along the a-axis found in IDPs (IBradley et al. 



19831 ) and /zm-sized polyhedral olivine cr ystals with various morphologies found in matrix of 
primitive meteorite ( jNozawa et al.ll2009l ). Hereafter, we refer to these fine crystals as cosmic 
crystals. It is an important issue to reveal the formation environment of these cosmic crystals 
for understanding the early history of the solar system. 

To reproduce the cosmic cryst als, evaporation and cond ensation experiments have been 
performed by many authors so far. iMysen fc Kushirol (119881 ) carried out experiments in the 
system composed of Mg2Si04-Si02-H2 in the pressure ranging from 10~ 4 to 10 4 dyn cm -2 and 
in the temperature ranging from 1620 K to 1920 K to determine phase relations of the system. 
In these experiments, condensation of MgSiC>3 and SiC>2 phases too k place from the Si-r i ch va - 
por, which was produced by incongruent vaporization of enstatite. iTsuchiyama et al.l (119881 ) 
analyzed these condensates with the use of an analytical transmission electron microscope 
(ATEM) and a scanning electron microscope (SEM) to compare the conden sates produced in 
the experiment with enstatite crystals found in IDPs (IBradley et al.lll983l ). They concluded 
that the characteristic textures of clinoenstati te found in IDPs were n ot reproduced in the 
evaporation and condensation experiments by lMysen fc Kushirol (119881 ). 



Recently, iKobatake et al.l (120081 ) carried out evaporation and condensation experiments 
to investigate a relationship between the growth conditions and morphologies of cosmic crys- 
tals condensed from highly supersaturated vapors. They used a sphere with forsteritic com- 
position (Mg 2 Si04) as an evaporation source and succeeded in reproducing various morpholo- 
gies observed in /xm-sized cosmic olivine crystals. They found that types of the morphology 
depend on the condensation temperature T c . Namely, they found bulky-types at T c > 1270 K, 
platy-types at 970 < T c < 1270 K, and columnar-needl e-types a t 770 < T c < 1090 K under 
the total pressure of 10 3 to 10 4 dyn cm -2 . Furthermore. lYamadal ( 120091 ) carried out the same 
experiments except that the evaporation source has enstatitic composition (MgSiOs) and suc- 
ceeded in reproducing enstatite whiskers elongated along the a-axis at 700 < T c < 1150 K. 
They also reproduced enstatite crystals of platy-types at 1150 < T c < 1300 K. Crys- 
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tals produced in their experiments have morphologies very similar to the cosmic crystals 
( IBradley et al.lll983t iNozawa et al.l 120091 ) . The condensation temperatures of these cosmic 
crystal analogs produced in the experiments aforementioned are substantially lower than the 
temperature of ~ 1400 K, at which forsterite and enstatite c ondense in equil ibrium from 
the solar nebula gas with the total pressure of 10 3 dyn cm -2 (I Grossman! 1 19721 . e.g.,). The 
reproduction of the morphologies of cosmic crystals in the experiments suggests that the 
cosmic crystals were formed in highly supercooled vapors through nucleation and successive 
crystal growth; the condensation temperatures of the cosmic crystals are lower than those 
expected from the equilibrium theory by a few hundred degrees or more. To produce such 
supercooled silicate vapor, one would require significant evaporation of silicate dust followed 
by rapid cooling of its vapor. 

As a possible site of formation of the cosmic crystals, we notice a localized bow shock 
produced by a planetesimal revolving in a highly eccentric orbit in a gas of the solar n ebula; 
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enter the bow shock, they do not evaporate significantly because of their large sizes (~ mm) 
but melt then cool and solidify to form chondrules. However, small dust particles of /im 
in size will evaporate c ompletely behind the shoc k front if the shock is strong enough, and 



produce silicate vapor (IMiura &: Nakamotoll2005l ). The vapor cools rapidly with expansion 



in the vicinity of the planetesimal in orbit, resulting in a highly supercooled state. 

In this paper, we examine formation of cosmic crystals in the regions behind planetesimal 
bow shocks. The point of our discussion is whether the silicate vapor realizes the highly 
supercooling in which cosmic crystals of various morphologies are produced, or not. We 
give an overview of our model in §|2j Section [3] describes the process of dust evaporation 
by shock-wave heating to evaluate the evaporation fraction of silicate dust particles in a 
wide range of the physical parameters. In §4J we investigate expansion of the silicate vapor 
behind the bow shock and estimate the cooling rate. Section [5] examines the properties of 
the condensates such as their particle sizes and morphologies expected from our model and 
compare the results with those obtained by the evaporation and condensation experiments. 
We discuss a comprehensive scenario on the formations of chondrules, cosmic crystals, and 
other materials in chondrites in 3BJ 



2. Outline of the Model 

As stated in Jfl a possible site for formation of fine cosmic crystals is a localized bow 
shock region associated with a planetesimal orbiting in an eccentric orbit. There arises 
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relative velocity between a planetesimal and the nebul ar gas both orbiting around the Sun, 
where the eccentricity of the planetesimal orbit is large. IWeidenschilling et al.l (119981 ) showed 
that the Jup iter mean motion reson ances can excite planetesimal eccentricities up to e ~ 
0.3 or more. iNagasawa et al.l (120051 ) analyzed the orbital evolution of terrestrial planetary 
embryos including the effect of the sweeping Jupiter secular resonance combined with tidal 
drag during dissipation of the protoplanetary gas disk. They found that the eccentricities of 
planetary embryos with mass of O.OIMe are excited up to e ~ 0.6 or more at maximum and 
oscillate around the mean value of e ~ 0.3 — 0.4 if Jupiter has the eccentricity of ej = 0.05 
(the current value is ej = 0.0485), where Me is the Earth mass. The relative velocity between 
the eccentric planetesimal and the circularly orbiting nebular gas is estimated to be v p ~ 
y/ e 2 + i 2 t>K, where e and i are, respectively, eccentricity and inclination of a planetesimal 
orbit and vk is its Keplerian velocity. For % <C e, we obtain 



10.3 



0.6/ V3 AU 



-1/2 



km s , 



where a is semi-major axis of a planetesimal orbit. The superso nic velocities relative to 



the n ebula gas produce bow shocks in front of the planetesimals (IHoodl Il998l ; ICiesla et al. 
20041 ) . The relative velocity of v p = 10.3 km s _1 for e = 0.6 at a = 3 AU is large enough 



to cause melting of mm -sized silicate dust aggreg ates (llida et al.l 1200 ll ) and evaporation of 
//m-sized dust particles (IMiura fc Nakamotoll2005l ) in a region of the asteroid belt, where the 
gas density is r in ~ 3 x 10 13 cm" 3 a t the midplane of the gas disk in the minimum mass solar 
nebula model (IHayashi et al.1119851 ). 



An outline of the model for cosmic crystal formation is illustrated in Fig. [TJ When 
the shock front formed by supersonic orbital motion of a planetesimal passes through a 
region in the nebula, the nebular gas is abruptly accelerated whereas /zm-sized dust particles 
tend to keep their initial position because of their relatively large inertia. As a result, 
the dust particles find that they are exposed to a high-velocity gas flow suddenly and are 
heated to their evaporation temperature if the relative velocity and the gas density are 
large enough. Evaporation of the /zm-sized dust particles in the post-shock region will be 
discussed in §3J The silicate vapor thus produced expands outward because its pressure is 
higher than that of the ambient unshocked region. The cooling associated with the expansion 
will produce supercooled silicate vapor. We discuss the cooling process of the silicate vapor 
and its cooling rate in §U The present model supposes that the cosmic crystals observed in 
meteorites are condensation products in the expanding silicate vapor supercooled behind the 
planetesimal bow shock. Their sizes, morphologies, and condensation temperatures strongly 
depend on the density of silicate vapor and the cooling rate. We shall show that various 
kinds of cosmic crystals observed in meteorites are formed in the cooling of the silicate 
vapor produced by planetesimal bow shocks (§5]). We shall point out that the present model 
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leads to simultaneous formation of chondrules and fine cosmic crystals ($6}, and that their 
formation is an inevitable consequence of formation of planetary systems. 



3. Evaporation of pm-Sized Dust Particles 
3.1. Evaporation fraction 

We carry out numerical simula tions of shock- wave heatin g by using a one- dimensional 



plane-parallel model developed by iMiura &: Nakamotol (120061 ). Actually, the structure of 



nebular gas around a planetesimal is not of one-dimen sional plane-paralle l structure. How- 



ever, the two-dimensional hydrodynamic simulation by iCiesla et al.l ( 12004 ) showed that the 
one-dimensional plane-parallel approximation was valid in the vicinity of a planetesimal we 
are concerned with, say, a few times of planetesimal radius i? p . We set a computational do- 
main along the x-axis to be — R p < x < R p , where the x-axis is parallel to the gas flow and 
x = at the shock front. I n this region, the shocked gas structure can be regarded as one- 



dimensional plane-parallel ( IHood et al.ll2009[ ). The simulations were carried out with varying 
the following input parameters: the planetesimal radius R p , the pre-shock gas number den- 
sity no, the dust-to-gas mass ratio £, and the shock velocity v s . In the simulations we set 
the ranges of parameters to be 1 < R p < 1000 km, 10 13 < n < 10 15 cm -3 , 0.01 < £ < 0.1, 
and 5 < v s < 60 km, respectively. We take particle radius to be = 1 /im as a typical size 
of fine dust particles; a scaling to other sizes is easily done with the use of Eq. fl3]). The case 
of £ = 0.1 is investigated to see the dependence of the evaporation fraction on £, although it 
will require significant settling or concentration of dust particles. 

Figure |2] shows the result for R p = 100 km, no = 10 15 cm -3 , v s = 8 km s -1 , and £ = 0.01. 
Panel (a) shows temperature profiles of the gas (T g , solid line), the dust particles (T d , dashed), 
and the radiation field (T ra( j, dotted) in the vicinity of the shock front (—1.0 < x < 1.0 km). 
Dust temperature increases rapidly just behind the shock front by gas frictional heating due 
to the velocity difference between gas and dust. In this stage, which we call the first stage, 
the dust temperature is determined by a balance among frictional heating, radiative cooling, 
and interaction with the ambient radiation field. The first stage ceases in a short period of 
time (less than 0.1 s in this case) because dust particles come to stop relative to the ambient 
gas. Panel (b) shows density profiles of dust (pa, dashed) and silicate vapor (p v , solid), which 
is produced by evaporation of the dust. It is found that the dust density pd increases behind 
the shock front (0 < x < 0.3 km) because of deceleration by the gas friction. On the other 
hand, the vapor density p v remains much smaller than pd, indicating that evaporation of the 
dust particles during the first stage is negligible because of the very short duration of the 
first stage. 
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Panel (c) shows temperature profiles over a wide region around the planetesimal. The 
relative velocity between the gas and the dust particles is almost zero in almost all the 
region shown here (x > 0.5 km), so the frictional heating does not work. However, the dust 
temperature is kept above 1500 K because of efficient collisional heating by the ambient 
hot gas; we call this stage of collisional heating the second stage. During the second stage, 
dust particles continue to evaporate gradually as is seen from the density profiles shown in 
panel (d). One should note that evaporation of dust occurs mainly in this stage. At the 
edge of the calculation zone (x = 100 km), the vapor density is p v = 1.76 x 10~ 10 g cm -3 , 
while the density of survived solid dust particles is pa = 1.73 x 10~ n g cm -3 . Therefore, the 
evaporation fraction rj defined by 

»? = -?- ( 2 ) 

Pd + Pv 

equals 0.91 in the case shown in Fig. [2j implying that 90 % of the dust mass evaporate during 
the second stage. 

Figure [3] summarizes the evaporation fraction 77 for various sets of values of the input 
parameters. Here, rj is shown as a function of dust temperature at the second stage. We take 
the dust temperature T^2 at the time when f re i/Vr = 0.1, where vt is root mean square of 
thermal velocities of the gas molecules and v Te \ is the velocity of dust particles relative to the 
gas; the result changes little even if we take v re i/vx = 0.05. Open circles in panel (a) show 
numerical results for all sets of the input parameters. Figure [3] indicates that rj increases 
rapidly with increasing The temperature dependence of rj is given by 

( Aa\ 3 

77 = 1-1 , 3 

V a d J 

where Aa is size decrease of a dust particle after finishing substantial evaporation and is 
given by 

a _ jcva,p(T d2 ) T d2 T d2 , . 

a ~ Pc H ~{dT/dt) Td2 1 } 

for a spherical dust particle (see Appendix |A] for the derivation). Here, p c is material density 
of the dust particle, j evap is the evaporation rate, i.e. mass of vapor evaporating from unit 
surface area of the particle per unit time, and H is latent heat of evaporation divided by the 
gas constant. Note that Aa is independent of the original size ad and r\ is small for a large 
dust particle. The factor 

At cvap = j^At (5) 
indicates an effective duration of evaporation during the cooling timescale defined by 



(6) 
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at T = Td2- Note that the numerical results are reproduced well by Eqs. (j3J) and @. This 
implies that the evaporation fraction 77 mainly determined by dust temperature in the 
second stage in spite that there are many other factors (Rp, n , v s , and £) that may affect 
evaporation of the dust particles behind planetesimal bow shock. 

The duration of substantial evaporation At evap is proportional to the cooling timescale 
At of the hot gas. The expression of At is very complex in general because it depends on var- 
ious physical processes such as vibrational/rotational transitions of CO and H 2 molecules, 
thermal dissociation of H 2 molecules, Lyman-a emission, and so forth. For a gas of the solar 
abundance, however, the major cooling process is L yman-q emission for T > 10 4 K and 



thermal dissociation of H2 molecules for T > 3000 K (IMiura fc Nakamotoll2005l . see Fig. 7). 
The timescale of cooling due to Lyman-a emission is shorter than ~ 100 s. Below 3000 K, 
the gas cools within a timescale of ~ 100 s du e to vibrational/rotational transitions of CO 



and H2O molecules (IMiura fc Nakamotoll2005l . see Fig. 7). The cooling timescale of the hot 



gas does not depend on the number density of the gas significantly. In the present case, the 
cooling timescale of At ~ 100 s reproduces the numerical results well as is seen from Fig. [3j 

Finally, let us examine the dependences of the evaporation fraction 77 on the parameters 
other than the temperature. Panels (b), (c), and (d) examine the dependence of the evap- 
oration fraction rj on R p , no, and £, respectively. In panel (b), n for R p = 1, 10, 100, and 
1000 km are plotted by different symbols to see the dependence of rj on R p . There seems 
no clear systematic dependence of 77 on R p even if we vary R p by three orders of magnitude. 
Panel (c) examines the dependence on n , the number density of pre-shock gas. There seems 
to be a slight trend that 77 decreases with increasing no but the dependence is unremarkable 
compared with the scatter of the data for each value of no. Panel (d) examines the depen- 
dence on the gas-to-dust mass ratio ranging from £ = 0.01 to 0.1 but we found no systematic 
trend of 77 on £, neither, within the plausible range of £. 



3.2. Analytic estimation of the dust temperature 

We have shown that the evaporation fraction 77 is determined mainly by the dust tem- 
perature Td2 in the second stage. However, one needs to elaborate numerical simulations 
to calculate T^ 2 . Instead, we derived an approximate analytic expression (1B2p of T^ 2 in 
Appendix [B] by considering the energy balance of a dust particle in the second stage. The 
analytic formula of T^2 will also be useful for calculating the dust temperature and its evap- 
oration in a planetesimal bow shock in general. 

Figure H] compares T^ 2 given by Eq. (1B2[) with that obtained from the numerical results. 
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It is found that both agrees with the difference less than ±50 K for T&2 < 1500 K. For 
^d2 ^ 1500 K, the numerical values of T d2 are systematically lower than those given by 
Eq. flB2p . The reason of the deviation is that the analytic estimation ignores decrease in 
the optical depth due to dust evaporation in the shocked region. Actually, the decrease in 
the optical depth weakens the intensity of the ambient radiation field, which heat the dust. 
In consequence, the dust temperature decreases and its evaporation is suppressed. This 
negative feedback taken into account in the numerical simulation results in the numerical 
value of Td2 lower than that of the analytic estimation. The deviation at > 1500 K, 
however, does not influence the estimation of the evaporation fraction rj much because r] ~ 1 
in any case at these temperatures as seen from Fig. [31 

Figure [5] shows the evaporation fraction r\ as a function of dust temperature T d2 as does 
Fig. El but in the horizontal axis is replaced by the one calculated by using Eq. (1B2[) . 
Although the scatter of the data plotted is larger than in Fig. [3l we see that the analytic 
formulae still reproduce the evaporation fraction 77. 



4. Expansion and Cooling of the Shocked Gas 

4.1. Equation of expansion 

When the hot gas in the shocked region cools down to the temperatures lower than 
~ 1500 K, dust particles re-condense from the vapor produced by evaporation of the original 
dust. In this subsection, we consider hydrodynamics and cooling of the expanding gas 
cloud to characterize the environment for formation of the cosmic crystals. Let us assume 
cylindrical expansion with initial radius Rq (see Fig. [I]). Initial radius of the shock ed region Rq 



is on the same order of magnitude as planetesimal radius i? p (jCiesla et al.ll2004j ). Neglecting 

(7) 



the expansion along the x-axis, the expansion velocity v r is described by 

dv r 1 dp 



dt p dr 

where p is the gas density and p is the gas pressure. We use a one-zone approximation and 
approximate v r and dp/dr as 

dR dp p 
*~dF' (8) 
where R is radius of the gas cloud at time t. We adopt a polytropic equation of state for the 
gas given by 

P ^)\ (9) 
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where po is initial gas density and 7 > 1 is a parameter relating to the polytrope index. The 
conservation of mass during the expansion is expressed as 

R 2 P = RW (10) 
Using Eqs. ([7]) to (ITUl) . we obtain the equation of expansion of the gas given by 

d 2 R 1 - 



dt 2 7 

with R = R/Rq and t = (Ro/c s o)t, where 



-fT 27+1 , (11) 



cso = J** (12) 
V Po 

is sound speed in the gas at t — 0. The dimensionless equations for expansion make it clear 
that the timescale of expansion of the shocked gas behind a planetesimal can be scaled by 
the sound- crossing time Rq/c s0 . 

Figure M shows the solutions of Eq. for the initial conditions of R = 1 and v r = 
(see Appendix [C]) . It is clearly seen that the expansion is separated into two phases; the 
acceleration phase, in which v r increases with time but R remains almost at the initial radius 
Rq, and the expansion phase, in which the shocked region begins to expand and v r almost 
equals a constant terminal velocity. The dashed curves in panel (b) show approximations of 
v r in the two limiting cases of t 1 and t — > 00 (see Appendix O) given by 

— 7- (c s0 t < Rq), 

7 t s0 . . 

I -[ 1/2 (I 3 ) 

(C s0 t > Rq), 



C S 



7(7 - 1) 



where t s0 = R /c s0 . 



4.2. Cooling rate of the shocked gas 

Using the relation 



T = ^(p/po) 7 " 1 = T R~ 2 ^ (14) 

given by Eqs. ([9]) and (flOl) . and T oc p/p, we obtain the time variation of the gas temperature 
T as 

-f-l"- 2 ( 1 ^) 1/2 ^ +1 [ 1 -^ 1, i 1/2 < 15 > 
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with the use of Eq. fICll) in Appendix O One sees from Eq. ffT5]) that the cooling rate 
—dT/dt as a function of R increases with increasing R at first, reaches a peak, and decreases 
in proportion to R~ 2l+1 . Figure [7] shows the cooling rate —dT/dt as a function of T which 
decreases monotonically with time. 

To evaluate the cooling rate, we need to specify a value of the initial temperature To. 
The gas temperature just behind the shock front could be higher than 2000 K or more 
depending on the Mach number v r /c S Q. However, even if the temperatures of the gas exceeds 
2000 K, it cools rapidly by dissociation of hydrogen molecules and is kept around 2000 K 
owing to the energy balance between re-form ation of hydrogen molecules by three-body 



reaction and their dissociation (jlida et al.ll200ll ). We set Tq = 2000 K to estimate the cooling 



rate around the condensation temperatures. To consider condensation through nucleation, 
on the other hand, we should refer to the cooling rate —dT/dt when the vapor becomes 
supersaturated. Taking the equilibrium condensation temperatures of T e = 1 444 K for 



forsterite and T e = 1349 K for enstatite for the total pressure of 10 3 dyn cm 2 (j Grossman 



19721 ) as a measure of estimating the condensation temperature, we have T c /T = 0.65 — 0.75 
and \dT/dt\ Tc = (0.25 - 0.35)T /(R /c s0 ) for 7 = 7/5 and 5/3 (see Fig. [7]). We set R to be 
planetesimal radius R p in what follows. In consequence, the cooling rate is estimated to be: 

fdT\ , . T 

- — ~ (0.25 - 0.35) 



\dt J Te -Rp/c s o 

* 2000 ( ~* ( ( c *° \ ! 

VI km/ V2000 Ky V3.7 km s- 1 / V ; 

It should be pointed out that cooling of the shocked gas given by Eq. (flEI) can be used 
so far as the pressure of the shocked gas p is much larger than the ambient gas pressure 
Pamb- The shocked gas pressure before the expansion is p ~ 100 p am b for the shock velocity 
of an H2 gas of 10 km s _1 . The gas temperature at th at time is ~ 2000 K as a result of the 



balan ce between H2 dissociation and its re- formation ( jlida et al.ll2001t iMiura &: Nakamoto 



20051 ). The pressure and temperature decrease by subsequent cylindrical expansion. When 
the temperature drops to the typical condensation temperature of ~ 1000 K, the radius 
of the cylinder is 2.4 times the initial one for adiabatic expansion, and the shocked gas 
pressure also decreases to ~ 1/ 10 of that before expansion. However, the gas pressure is still 
higher than p am b by an order of magnitude. Therefore, Eq. (fl6|) is applicable throughout the 
expansion phase of interest including the time of condensation. 

We focus here the adiabatic expansion because the radiative losses are negligibly small 
for small shocks as is shown below. Main coolants of the nebula gas at 2000 K are vibra- 
tional emissions of CO and H2O molecules. The cooling timescale due to these vibrational 
emissions was estimated to be ~ 100 sec, which does not significantly depend on the gas 
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density ( IMiura fc Nakamotol 120051 ). On the other hand, the cooling timescale due to the 
adiabatic expansion behind a planetesimal is shorter than ~ 100 sec for planetesimal radius 
of < 100 km (see Eq. ( TTBj) ). Therefore, the shocked gas cools by the expansion before the 
vibrational emissions remove the thermal energy significantly. The radiative losses might 
work for large shocks (> 100 km) because the expansion takes longer time. However, a large 
optical depth for these emissions resulting from large shocks prevents the radiative losses 
from being efficient. 



4.3. Possibility of chondrule formation 



In the formation of chondrules, their cooling rate during solidification is one of the key 
physical quantities. According to the planetesimal bow shock model, the cooling rate was 
estimated to be ~ 10 3 K hr" 1 for planetesimal radius R p = 10 00 km, > 10 4 K hr -1 for 
R p = 100 km, and > 10 5 K hr" 1 for R p = 10 km Jrlood et al.ll2005h . On the other hand, the 
cooling rate of the shocked gas calculated from Eq. (1161) is 7 x 10 3 , 7 x 10 4 , and 7 x 10 5 K hr -1 
for R p = 1000, 100, and 10 km, respectively. Although Eq. ffl6|) is not a cooling rate of a 
chondrule itself but of the shocked gas strictly speaking, we note that both estimations of 
the cooling rates are co mparable; this is because the cooling of chondrules is regulated by 
that of the shocked gas (llida et al.ll200ll ). Therefore, Eq. ffTB]) measures the cooling rate of 
chondrules. 

A widel y accepted range o f the cooling rate of chondrules at solidification is ~ 10 — 
1000 K hr -1 ( jHewins et al.ll2005l . and references therein), which is much slower than that pre- 
dicted by Eq. ( TTB|) . However, we consider that this disagreement does not necessarily exclude 
planetesimal bow shocks as a chondrule formation site. In fact , some authors asser t rapid 
cooling rates, which are in the range estimated from Eq. f JT6|) . lYurimoto &: Wassonl (120021 ) 
proposed that rapid cooling (~ 10 5 — 10 6 K hr -1 ) was necessary to accoun t for the observed 



Fe-Mg and O-isotopic exchange in a CO-chondrite type-II chondrule. iWasson &: Rubin 
( 120031 ) proposed that very thin overgrowths on some relict grains in chondrules must have 
been formed by the rapid cooling. The crystallization experiments of a melt droplet by a levi- 
tation method succeeded in repr oducing chondrule-solidification textures in the e xperimental 
conditions of the rapid cooling (ITsukamoto et al.lll999t iNagashima et al.l 120061 ) . Although 
the rapid co oling scenario does not seem to have been accepted widely to the meteoritic 
community ( jHewins et al.l 120051 ). there has been no definite evidence that rejects the rapid 
cooling scenario completely. Therefore, we consider that the planetesimal bow shock is still 
one of the possible models to be studied for chondrule formation. 
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5. Formation of Cosmic Crystals 

5.1. Cooling parameter A for homogeneous nucleation 

Cosmic crystals condense in the course of cooling of the vapor produced by a planetes- 
imal bow shock. When almost all dust particles evaporate by the bow shock, there is no 
solid surface available on which the supersaturated vapor condenses. In this case, cosmic 
crystals are formed through homogeneous (spontaneous) nucleation. In homogeneous nucle- 
ation, condensation does not begin when the cooling vapor becomes saturated but begins 
effectively after the vapor becomes supersaturated to a certain degree. 



Yamamoto fc Hasegawal (119771 ) and lDraine fc Salpeten (119771 ) formulated a grain forma- 



tion process though homogeneous nucleation in a vapor and derived analytical expressions 
of a typical size of grains and their actual condensation temperature T c in a supercooling 
state as functions of two dimensionless parameters. One is a cooling parameter defined by 

where v co \\ is collision frequency of vapor molecules in thermal motion, tx = T e / (—dT/dt) T ^ 
is cooling timescale of a vapor at T = T e with T e being equilibrium condensation temperature, 
and H is latent heat of condensation divided by the gas constant and equivalent to that of 
evaporation (see Appendix [A} . Note that T e (> T c ) is a temperature at which a vapor and 
a bulk condensate co-exit in chemical equilibrium and approaches T c as tx gets so long that 
the equilibrium between the vapor and the condensate is realized. Grain size aa is mainly 
determined by A an d is roughly given by an /an ~ 0.1 A for A ^> 1, where a is the radius 



of a vapor molecule ( lYamamoto fc Hasegawal ll 9 771 ). In Eq. (fTTj) . v co \\ is calculated from the 



vapor density p v , and fa from the cooling rate of the vapor (see Appendix [D]). The other 
parameter is a dimensionless surface tension defined by 

r = T^, (18) 

where 7 S is surface tension of a condensate and an = (3/i c m a /47rp c ) 1 / 3 (i.e. equivalent radius 
of a sphere whose volume equals the volume of a unit cell of the condensate) with fi c being 
molecular weight of a unit cell of a condensate, m a = 1.66 x 10~ 24 g is atomic mass unit, 
and p c its bulk density. A degree of supercooling AT = T e — T c i s mainly determined by 



the p arameter T and is approximately related to T as AT oc T 3//2 (lYamamoto fc Hasegawa 



1973). 



Figure [8] shows a relation between A and the evaporation fraction rj. Each of the plots 
indicates rj calculated in £J3] for a given set of values of the parameters, while A is calculated 
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from Eq. (ID2j) . All panels indicate the trend that A increases with 77. This is simply because 
the larger degree of evaporation of pre-existing dust is, the larger amount of the vapor is 
produced, which in consequence provides favorable conditions for homogeneous condensation 
of cosmic crystals. Note that homogeneous condensation is possible only if A > 1; otherwise, 
the vapor is too tenuous for condensation to occur. Panel (a) shows the results of the 
calculations for all of the parameter sets, indicating that there appear many cases of A > 1 
for 77 > 10 -4 . Even the cases of A as large as 10 5 are realized for complete evaporation 
(77 ~ 1) of pre-existing dust. The contribution to the vapor production comes mainly from 
/jm-sized dust particles if their size distribution is steeper than a^ 2 . The presence of many 
cases of A > 1 implies that condensation of cosmic crystals through homogeneous nucleation 
behind planetesimal bow shocks is possible for 77 > 10 -4 . We note that a variety in the 
A-values suggests formation of various kinds of cosmic crystals. The panel (b) displays the 
dependence of 77 and A on the planetesimal radius R v . From panel (b), one sees that the 
homogeneous condensation occurs hardly except for 77 ~ 1 for a bow shock produced by small 
planetesimals of R p = 1 km but occurs almost always for a planetesimal of R p = 1000 km 
even if the evaporation is not so significant (77 > 10" 4 ). 



5.2. Size and morphology of cosmic crystals 

Figure M displays typical size of condensed particles and supercooling AT in terms of 
A and V. The supercooling AT in the vertical axis is normalized by the equilibrium condensa- 
tion temperature T e . Each solid curve shows the relation between a y and AT for a constant 



value of T, and dashed lines combine points for the same value in A (jYamamoto fc Hasegawa 



19771 ). The grayed region indicates a parameter range expected from the planetesimal bow 
shock. The possible range of A was discussed in §5.11 The values of T, on the other hand, 
have uncertainties because of short of the experimental data for the surface tension 7 S of 
forsterite and enstatite. For forsterite, 7 S is me asured to be 1280 erg cm -2 in vacuum for a 



{010} surface and larger values for other ones (Ide Leeuw et al.l 120001 ). which corresponds to 
T ~ 30 or more. For enstatite, there are no reliable data of surface tension. We assume the 
similar value as that of forsterite. In the calculations, we take 10 < T < 60 for safety. It 
should be noted that the sizes Oqo and the super coolings AT revealed from the analyses of 
a variety of cosmic crystals are included in the region realized by planetesimal bow shocks. 
Let us discuss in more detail the formation conditions of each of the cosmic crystals shown 
in Fig. El 
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5.2.1. Enstatite whisker and platelet 



The experiment by lYamadal (120091 ) showed that formation of enstatite whiskers elon- 
gated toward the a-axis required the degree of supercooling of 0.15 < AT/T C < 0.48. They 
also reproduced platy-type enstatite crystals at 0.04 < AT/T C < 0.15. It is interesting to 
note that the whisker has larger T than the platy-type, although precise values of their sur- 
face tension are unknown. Typical size of the enst atite crystals is ~ 0. 1 - 1 /im, which size 
is similar to that of natural samples found in IDPs ( jBradley et al.lll983l ). A set of these con- 
ditions is shown by the red region in Fig. [9j indicating that enstatite whiskers and platelets 
can be formed by planetesimal bow shocks of 10 3 < A < 10 4 . This range of A is realized if 
the bow shocks are produced by planetesimals of intermediate size (_R P ~ 100 km) and lead 
to almost complete evaporation of the original dust (77 ~ 1). If the amount of the silicate 
vapor is small leaving a large amount of dust particles that survived evaporation < 1), on 
the other hand, the vapor will condense onto the dust surface. This case yields other types 
of thermally-processed particles observed in chondritic meteorites (see 



5.2.2. Olivine crystals with various morphologies 



Kobatake et al.1 (120081 ) examined supercooling AT required for formation of olivine crys- 



tals by a laboratory experiment. They showed that bulky-type olivine crystals were repro- 
duced at AT/T C < 0.12, the platy-type at 0.12 < AT/T e < 0.33, and the columnar-needle- 
type at 0.24 < AT/T C < 0.47. As was so for enstatite, the needle-type has larger Y than the 
platy-type; the bulky type has the smallest Y. Typical size of the condensates is ~ /an, 
which size is close to those of the natural samples found in the matrix of Allende meteorite 
( iNozawa et al.ll2009l ). The green region shows the supercooling AT/T e and the sizes o m for 
these fine olivine crystals, indicating that these particles can be formed by planetesimal bow 
shock of 10 4 < A < 10 5 . This condition is realized by the bow shocks produced by relatively 
large planetesimals (R p ~ 1000 km) associated with almost complete evaporation of the 
original dust particles. 



5.2.3. Ultra- fine particles 



Toriumil (119891 ) observed fine particles in the matrix of Allende meteorite using a SEM 



and a TEM and measured their sizes a^. The observed size distribution could be reproduced 
by a log-normal one for 1 < ad < 10 nm with its peak at a = 5 nm and by a power law for 
aa > 10 nm. We display the size range of the ultra- fine particles by the blue region in Fig. |9j 
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The size range suggests that A ~ 10 — 100 is a plausi ble condi t ion fo r formation of ultra-fine 
particles. This is in agreement with the conclusion of iToriumil ( 119891 ) that ultra-fine particles 
seem to have been formed by condensation from a vapor far from equilibrium in the early 
solar nebula. The present model implies that ultra-fine particles were formed by bow shocks 



produced by much smaller planetesimals (_R P 



10 km) than those producing /xm-sized 



cosmic crystals, associated with almost complete evaporation. The formation condition of 
A ~ 10 — 100 also realizes for large planetesimals (R p > 100 km) and small evaporation 
fraction (77 ~ 10~ 4 — 10~ 3 ), however, in this case the ultra-fine particles generated from the 
vapor are very rare because of the tiny evaporation fraction. 



5.3. Heterogeneous condensation for incomplete evaporation 

We discussed formation of cosmic crystals through homogeneous nucleation in §5.11 and 
§5. 21 assuming that almost all dust particles evaporate by a planetesimal bow shock. There is 
an opposite case that a substantial fraction of the dust particles survives against evaporation 
and acts as seed nuclei and that condensation occurs through nucleation on their surfaces 
(heterogeneous condensation). Which type of condensation actually occurs depends on the 
total surface area of dust particles available for heterogeneous nucleation. We shall show 
below that both types of condensations can occur depending on the radii of planetesimals 
generating bow shocks and on the evaporation fraction. 

In homogeneous nucleation, condensation does not begin when the cooling vapor be- 
come saturated but begins effectively after the vapor becomes supersaturated to a certain 
degree. Namely, there arises s ome induction time t inc i after the vapor becomes saturated 



(jYamamoto fc Hasegawal Il977f ). The induction time is related to the cooling timescale 



X J f n no a n\ ( 



t T = T e /(-dT/dt) Tc as 

^"i^ 4 ^' 08 - 4 - '^ s - (19) 

where xj = 2 — 70 for situations we consider in this paper (A = 1 — 10 5 and T = 10 — 60 as 
explained in §5.11 and §5.2p . The time inter vals required for nucleation and growth is about 



ten times shorter than the induction time (jYamamoto fc Hasegawal 119771 ). Therefore, ti n d 



represents a typical timescale for dust formation through homogeneous nucleation after the 
vapor becomes saturated. 

In heterogeneous nucleation, on the other hand, we estimate its timescale by using the 
adhesion timescale, which provides an underestimate of the timescale of heterogeneous con- 
densation because it ignores the induction time for heterogeneous nucleation. The adhesion 



- 16 - 



timescale t a d is the one during which most of the vapor molecules sticks onto the surface 
of dust particles. For silicate condensation, we regard SiO molecule as a key species that 
controls the rate of condensation (see also Appendix [D]). Denoting the radius of the dust 
particles by a^, the adhesion timescale is estimated to be: 

= (g^spy* g 5Q 1(1 _ \ / 10-" g cm- N 

3a s p d V k B T e J VP 111 / V Pd + Pv / 

where pa = (1 — r?)(p d + Pv) is density of dust particles surviving in the post-shock region 
against evaporation, a s is sticking probability of vapor molecules onto the dust surface, and 
PsiO — 44 is molecular weight of SiO. One should note that, in Eq. (1201) . the factor 3pd/Pc Q d 
indicates total surface area of the dust particles per unit volume and (8fcB^e/ 7r PsiO m a) 1//2 is 
mean thermal velocity of SiO molecules. 

Homogeneous nucleation takes place if ti n d < £ a d- This condition is satisfied when the 
planetesimal radius R p is relatively small (R p < 10 — 500 km), or there are few survived dust 
particles because of significant evaporation (77 ~ 1). In this case, cosmic crystals condense 
directly from the vapor. In contrast, heterogeneous condensation becomes effective if the 
planetesimal radius R p is large (R p > 10 — 500 km) and a substantial fraction of dust 
particles survives against evaporation (i| < 1). We shall discuss generic relations between 
cosmic crystals and chondrules in § [0J in detail. 



6. Summary and Discussion 



Chondritic meteorites are composed of materials that have been experienced thermal 
processing of various degrees in the early solar nebula. These materials includ e chondrules, 
fine-grained rims on chondrules and interchondrule matrix ( lAlexanderl Il995l ) , and cosmic 
crystals discussed in the previous section. In this section, we discuss how the planetesimal 
bow shock scenario explains the formations of these chondritic materials. 

A planetesimal bow shock was originally proposed as a possible site for chondrule forma- 
tion (lHoodlll998l ). Ilida et al.l ( 1200 ll ) showed that millimeter-sized dust aggregates (chondrule 
precursors) are heated and melt behind a shock front if the shock velocity and the pre-shock 
gas density are in an appropriate range. Complete evaporation hardly occurs for chondrule 
precursors because of their large size (see Eq. <^fy). Their contribution to the vapor pro- 
duction is negligibly small compared with that of /zm-sized dust particles for the dust size 
distribution steeper than aj 2 . Large molten chondrule precursor dust survives against evap- 
oration, cools and solidifies to form chondrules. In contrast, /im-sized particles evaporate 
significantly in the hot gas behind the planetesimal bow shock and produce silicate vapor. 
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The silicate vapor cools rapidly behind the bow shock and becomes supersaturated, leading 
to condensation to produce various kinds of materials observed in chondritic meteorites and 
IDPs. 

The condensed materials exhibit a wide variety in morphologies and sizes depending on 
their formation conditions such as the cooling rate and the evaporation fraction of //m-sized 
dust particles. The cooling rate is inversely proportional to the size of a planetesimal that 
produces a bow shock (see Eq. (fTBl) ). thus decreases with time on average, namely, with 
growth of planetesimals. The evaporation fraction rj changes by many orders of m agnitude 



in the range of the shock conditions realized in early solar nebula (llida et al.ll200ll ). 



Figure [10] summarizes condensation products in the course of the planetesimal growth. 
At the early stage of 1 < R p < 10 km, the vapor produced by small planetesimals cools 
so rapidly that the cooling parameter is A < 10 3 , which is realized for 10 -2 < rj < 1 (see 
Fig. [8b). Condensation of the vapor through homogeneous nucleation for 1 < A < 10 3 leads 
to formation of nm-sized ultra-fine particles as observed in the matrix (see from Fig. [9]). 
Furthermore, the results of §5.31 indicate that the rapid cooling prevents heterogeneous con- 
densation on survived dust particles because of rapid consumption of the vapor by homoge- 
neous condensation to form ultra-fine particles. To summarize, most of the vapor condensed 
to the ultra-fine particles at the early stage of planetesimal growth. When planetesimals 
grow up to a several 100 km or more, condensation occurs through both homogeneous and 
heterogeneous nucleations. If almost all of the small dust particles evaporate (rj ~ 1) behind 
the bow shock, /zm-sized euhedral silicate crystals condense through homogeneous nucle- 
ation. The cooling parameter is 10 3 < A < 10 5 for r\ ~ 1 (see Fig. [8b) for planetesimals 
of 100 < R p < 1000 km. This situation leads to condensation of enstatite whisker elon- 



gated to a-axis as found in IDPs ([Bradley et all Il983l) and polyhedral olivine crystals as 



found in the matrix of Allende meteorite ( jNozawa et al.ll2009l ) (see Fig. [H]). Bare chondrules 



without fine-grained rims could also be formed in this case. On the other hand, if many of 
the dust particles survive against evaporation (rj <C 1) and suffer partial evaporation, the 
vapor condenses heterogeneously onto the survived dust particles, resulting in the formation 
of other kinds of meteoritic materials. The vapor condensed heterogeneously on chondrules 
already solidified could form fine-grained rim on their surfaces. The survived /zm-sized dust 
particles would also be covered with materials condensed from vapor, and would accumulate 
as fine-grained interchondrule matrix in chondritic meteorites after that. It is worth noting 
that the partial evaporation of dust particles would lead to elemental fractionation. The 
fractionated vapor rich in volatile elements re-condensed within a short period of time (see 
Eq. ( I20p ) on the survived dust particles. A fine-grained rim of a chondrule and a fine-grained 
interchondrule matrix thus produced would have elemental composition complementary to 
that of the chondrule. The composition of the whole particle should be the same as that of 
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the original dust particles before evaporation according to the present model. This is consis- 
tent with the genetic relationship among chondrules, interchondrule matrix, and fine-grained 
rims that these components either formed from a common source material, are pr oducts of 



the same process, or have exchanged materials during formation (jHuss et al.ll2005l ) 



In summary, the planetesimal bow shock model can provide a comprehensive scenario 
for the formation of various cosmic crystals and other materials observed in chondritic mete- 
orites. Because the heating events happened in a localized region of the shocked gas within a 
short period of time, one may expect that a series of the thermal processing, heating, evapo- 
ration, and condensation completed in a closed-system. The scenario is in harmony with the 
genetic relationship suggested by the analyses of chondritic meteorites and IDPs that these 
are produced in the course of the processing from a common source material together with 
exchanges of the materials during their formation. 
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A. Size Decrease of a Particle by Evaporation and the Evaporation Fraction 

Let us consider evaporation of a spherical dust particle of initial radius ad- The evapo- 
ration fraction 77 is given by 

3 T d 9 , ( Aa\ 3 , k . 

77 = r / Air a 2 da =1-1 , (Al) 

47ra3y ad „ Aa \ a d J 

where Aa is decrease in radius due to evaporation. Equation (1A1I) indicates clearly that 
77 depends only on Aa/a^, the ratio of the size decrease to the initial size. Since the size 
decrease Aa due to evaporation is independent of the particle radius a d except through a 
slight dependence of the dust temperature on a d , the increase in the particle radius simply 
causes the decrease in 77 according to Eq. flAlj) . For example, even if /im-sized dust particles 
evaporate almost completely (77 = 0.999), chondrule-sized particles (od = 500 /im) evaporate 
by only a small fraction of rj ~ 5 x 10~ 3 . We carry out the calculations for ad = 1 /im, but 
one can evaluate 77 for other a d by using Eq. flAlj) . 
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With the use of the evaporation rate jevap, the size decrease Aa by evaporation during 
cooling from temperature T x to Tf is expressed by 



Aa = — I j 

Pc 



cvap 



dt= — 



Jevap 

(T) 

p c 7t; dT/dt 

The evaporation rate as a function of temperature T behaves as 

jcvap(^) = const ■ T p exp(-H/T), 



dT. 



according to the Hertz-Knudsen equation (INagahara et al.l Il996l ; iMiura et al.l 12002 



(A2) 



(A3) 



see 



where (3 ~ —1/2 is a constant. Here, H = L cvap /R g3bS is latent heat of evaporation in 
units of temperature, where L evap is that in units J mol -1 and i? gas is the gas constant in 
units of J K _1 mol -1 . In the present case, we are concerned with evaporation of forsterite 
(Mg 2 Si0 4 ), for which L evap = 1.58xl0 13 J mol -1 . This leads H = L eva , p /6R gas = 3.17xl0 4 K, 
where the factor of 6 results from the stoichiometric coefficients of the chemical reactions at 



evaporation ( IMiura et al.ll2002l see Eq. (38)). According to the measurement of evaporation 
rate of forsterite, j e ^, p de pends also on the partial pressure pn 2 of ambient hydrogen molecule 
( ITsuchiyama et al.lll998l ). but we may take the pressure at T = 7] in Eq. flA2j) because Aa 
is determined by the physical conditions at T = T ; as will be seen below. Integration on 
the RHS of Eq. 0A2p can be performed by noting that e~ H l T is a rapidly varying function 
compared to the remaining function in the integrand. Integrating by part and remaining the 
term of order T\jH <^ 1, one obtains 



Aa = Jevap(Ti) ^At, 
Pc H 

where the contribution from the upper limit of the integral is negligible. Here 



(A4) 



At 



-dT/dt) Ti 



(A5) 



is cooling timescale of dust particles at T = TJ, for which we take dust temperature in 
the second stage. 



B. Dust Temperature behind Shock Front 

We give here an analytic expression that gives in good approximation of the dust tem- 
perature in the post-shock region after the relative velocity between the gas and dust particles 
is almost damped (the second stage). The dust temperature at this stage, Td2, is determined 
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by the energy balance between collisional heating by the ambient hot gas and the radiative 
cooling: 



1 7 + 1 fjk B T 



47 — 1 \7Tp,m s , 



f\ 1/2 



n'k B (T' - T d2 ) + a SB (T r 4 ad - T d 4 2 ) = 0, 



(Bl) 



where T' is post-shock gas temperature, n' is post-shock gas number density, T rad is ambi- 
ent radiation temperature, \i is mean molecular weight of the gas, and <tsb is the Stafan- 
Boltzman n constant. Here, we appr oximated the emission and absorption coefficients to 
be unity (IMiyake fc Nakagawal Il993l ). Since T ra d ~ T d2 in the second stage as seen from 
Fig. [2(c), one obtains 



T, 



(12 



^kvT'/ixpm^n'kvT'/A + 4^ 



rad 



(B2) 



7r {jk B T'/7Tfim a ) VVAfe/4 + 4as B T r 3 ad 

from Eq. ( IBip by using the approximation that (T rad — T d2 ) ~ 4T r 3 ad (T rad — T d2 ), where 
7r = ( 7 + l)/(7-l). 

In Eq. ( 1B2j) . the post-shock gas number density n' is given from the Rankine-Hugoniot 
relation and the almost isobaric condition for the post-sho ck gas in the one-dimensional 
plane-parallel geometry ( jSusa et al.lll998l ; iMiura et al.ll2002l ) by 

2 



n 



(B3) 



7 + I^bT'' 

where po = pm^no is the gas density in the pre-shock region and v s is the shock velocity. 

The post-shock gas temperature T' in Eq. ( IB2[) may be obtained by using the Rankine- 
Hugoniot relation. However, we have to pay attention that, at high temperatur es of T' > 



2000 K, the gas cools very rapidly due to the dissociation of hydrogen molecules ( llida et al. 
200ll ) . Therefore, we set T' as 

2(7-1) pm & vl 



T' 



mm 



(7 + 1) 2 h 



2000 K 



(B4) 



The radi ation temperature at the shock front is given by taking the blanket effect into 
account ( IMiura fc Nakamotoll2006l ) as 



rad 



3r, 



4a SB 2 



pre / 3 

Pov, 



1/4 



(B5) 



where r pre is optical depth of the pre-shock region and / is the fraction of the gas energy flux 
that returns upstream in the form of radiation. We set / = 0.5 for simplicity. The optical 
depth r pre is estimated to be 

3£po r 



pre 



4a d p c 



(B6) 
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where £ is dust-to-gas mass ratio in the pre-shock region and L s is dimension of the pre-shock 
region, in which the dust particles contribute to the blanket effect around the shock frontQ. 
We set L s = R p in this study (see §3.ip . 

Equation ( 1B2|) together with Eqs. ( 1B3j) . (1B4j) . and (1B5j) is an analytic expression of the 
dust temperature in the second stage. It should be noted that these equations include all of 
the input parameters L s = R p , n = p / ^m a , £, and v s . 



C. Solutions of the Equations of Expansion 



Integrating Eq. ( ITT]) from t — to t after multiplying dR/dt = v r on both sides, we 
obtain the expansion velocity to be: 



dR 
dt 



1 _ £-2(7-1) 
7(7 - 1) 



1/2 



(CI) 



for the initial conditions of v r = and R = 1. Equation ( IClj) is integrated further to yield 
the radius R as a function of time t as 



t 



^7(7 - 1) 



R 



dy 



^1-2,-2(7-1)- 



(C2) 



The right-hand side of Eq. (1C2j) may be expressed by the hypergeometric function but 
numerical integration is more practical to get the results, which are shown in Fig. [6](b) by 
solid curves. 



In the limits of £ <C 1 and t — > 00, the velocity is approximated to be: 

(t«l), 



t 

1 



7(7 - 1) 



1/2 



(C3) 



(t -)• 00) 



Both approximations are shown in Fig. [6]^b) by the dashed lines. The time at the intersection 
at which the two limiting approximations cross each other, is given by 



7 



7-1 



1/2 



(C4) 



1 L S corresponds to x m in Miura & Nakamoto (2006). 
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At the intersection, the radius and expansion velocity are given by 



R* 



2 7 - 1 
2(7-1)' 



1/2 



.7(7 - 1). 



(C5) 



D. Evaluation of cooling parameter A for multi-component evaporation 

In this paper, we are concerned with evaporation and co ndensation of forsterit e (A/h^SiC^), 
in which Mg and SiO should be considered as vapor species (INagahara et al.lll996l ). For deal- 
ing with nucleation of a multi-component system, we adopt the key species approximation 
that the rates of nucleation and grain growth are controlled by one chemical species (key 
species) that h as the least collision frequen cy among the major vapor species that condense 
into the grain ( Kozasa fc Hasegawal 1987). The condi tions for the key-species approxima- 
tion to hold were examined by lYamamoto et al.l ( 1200 ll ) in formulating theory of nucleation 
involving chemical reactions. 

In Eq. (I17p , the collision frequency of vapor molecules of mass /i v ^ a and number density 
n v is given by v co \\ = nala s n v ^8k B T e / 7r/i v m a , where q s is sticking probab i lity, u v is mean 
molecular weight of the vapor molecules. Following Kozasa fc Hasegawal ( 119871 ). we take 
SiO molecules as the key species of silicate condensation. This implies that n v = nsio and 
/-tv = HsiO — 44. Using Eqs. (TTBT) and (TT2"j) . we obtain 



A 



4ira 2 a s n SiO R p f /1 
H/T c -1 \fisio 



1/2 



To 



3/2 



1 

7T7 



1/2 



(0.25 -0.35)' 



(Dl) 



where fi is mean molecular weight of the gas. The value of H/T e is estimated to be H/T c — 1 ~ 
20 for H ~ 3 x 10 4 K for forsterite; enstatite yields the similar value. The value of ao is given 
by a = 2.6 A. In consequence, A is evaluated roughly to be: 



A ~ 400 



Rr 



100 km 



10 10 g cm 3 



(D2) 



for a s = 1. 



REFERENCES 

Alexander, C. M. O'D. 1995, Geochim. Cosmochim. Acta, 59, 3247 
Bradley, J. P., Brownlee, D. E., & Veblen, D. R. 1983, Nature, 301, 473 



-23 - 



Ciesla, F. J., Hood, L. L., Weidenschilling, S. J. 2004, Meteorit. Planet. Sci., 39, 1809 

de Leeuw, N. H., Parker, S. C, Catlow, C. R. A., & Price G. D. 2000, Phys. Chem. Minerals, 
27, 332 

Draine, B. T., & Salpeter, E. E. 1977, J. Chem. Phys., 67, 2230 
Grossman, L. 1972, Geochim. Cosmochim. Acta, 36, 597 

Hayashi, C. K., Nakazawa, K., & Nakagawa, Y. 1985, in Formation of the solar system, ed. 
D. C. Black k M. S. Matthews, M.S. (Univ. of Arizona Press, Tucson), 1100 

Hewins, R. H., Connolly, Jr., H. C, Lofgren, G. E., & Libourel, G. 2005, in Chondrites 
and the Protoplanetary Disk, ed. A. N. Krot, E. R. D. Scott, & B. Reipurth (San 
Francisco: Astronomical Society of the Pacific), 286 

Hood, L. L. 1998, Meteorit. Planet. Sci., 33, 97 

Hood, L. L., Ciesla, F. J., & Weidenschilling, S. J. 2005, in Chondrites and the Protoplane- 
tary Disk, ed. A. N. Krot, E. R. D. Scott, & B. Reipurth (San Francisco: Astronomical 
Society of the Pacific), 873 

Hood, L. L., Ciesla, F. J., Artemieva, N. A., Marzari, F., & Weidenschilling, S. J. 2009, 
Meteorit. Planet. Sci., 44, 327 

Huss, G. R., Alexander, C. M. O'D., Palme, H., Bland, P. A., & Wasson, J. T. 2005, in 
Chondrites and Protoplanetary Disk, ed. A. N. Krot, E. R. D. Scott, & B. Reipurth 
(San Francisco: Astronomical Society of the Pacific), 701 

Iida, A., Nakamoto, T., Susa, H., & Nakagawa, Y. 2001, Icarus, 153, 430 

Kobatake, H., Tsukamoto, K., Nozawa, J., Nagashima, K., Satoh, H., & Dold, P. 2008, 
Icarus, 198, 208 

Kozasa, T., & Hasegawa, H. 1987, Prog. Theor. Phys., 77, 1402 

Miura, H., Nakamoto, T., & Susa, H. 2002, Icarus, 160, 258 

Miura, H., & Nakamoto, T. 2005, Icarus, 175, 289 

Miura, H., & Nakamoto, T. 2006, ApJ, 651, 1272 

Miyake, K., & Nakagawa, Y. 1993, Icarus, 106, 20 

Mysen, B. O., & Kushiro, I. 1988, American Mineralogist, 73, 1 



Nagahara, H. k Ozawa, K. 1996, Geochim. Cosmochim. Acta., 60, 1445 
Nagasawa, M., Lin, D. N. C, k Thommes, E. 2005, ApJ, 635, 578 

Nagashima, K., Tsukamoto, K., Satoh, H., Kobatake, H., k Dold, P. 2006, J. Crys. Growth, 
293, 193 

Nozawa, J., Tsukamoto, K., Kobatake, H., Yamada, J., Satoh, H., Nagashima, K., Miura, 
H., k Kimura, Y. 2009, Icarus, 204, 681 

Susa, H., Uehara, H., Nishi, R., k Yamada, M. 1998, Prog. Theor. Phys., 100, 63 

Toriumi, M. 1989, Earth Planet. Sci. Lett., 92, 265 

Tsuchiyama, A., Kushiro, I., Mysen, B. O., k Morimoto, N. 1988, Proc. NIPR Symp. 
Antarct. Meteorites, 1, 185 

Tsuchiyama, A., Takahashi, T., k Tachibana, S. 1998, Mineralogical Journal, 20, 113 

Tsukamoto, K., Satoh, H., Takamura, Y., k Kuribayashi, K. 1999, Antarct. Meteorites, 24, 
179 

Wasson, J. T. k Rubin, A. E. 2003, Geochim. Cosmochim. Acta, 67, 2239 

Weidenschilling, S. J., Marzari, F., k Hood, L. L. 1998, Science, 279, 681 

Yamada, J. 2009, Master's Theses, Department of Earth and Planetary Science, Tohoku 
University 

Yamamoto, T., k Hasegawa, H. 1977, Prog. Theor. Phys., 58, 816 
Yamamoto, T., Chigai, T., Watanabe, S., k Kozasa, T. 2001, A&A, 380, 373 
Yurimoto, H. k Wasson, J. T. 2002, Geochim. Cosmochim. Acta, 66, 4355 



This preprint was prepared with the A AS IATgX macros v5.2. 



-25 - 




Fig. 1. — Outline of the model for cosmic crystal formation. The nebular gas and precursor 
silicate dust come from the left side of the planetesimal and pass through the shock front 
produced by a planetesimal orbiting at supersonic velocity in the nebular gas. They are 
heated behind the shock front, and evaporation of the dust particles produce silicate vapor, 
which comes mainly from evaporation of /zm-sized particles. The grayed region indicates 
existence of the vapor. Pressure, density, and temperature of the gas in the shocked region 
just before the expansion (t = 0) are denoted by p , p , and T , respectively. The shocked 
region has higher pressure than the ambient region and expands vertically with velocity 
v r . R(t) is radius of the gas cloud at time t. The silicate vapor cools with expansion 
and becomes supercooled. Cosmic crystals condense from the cooled vapor after the vapor 
becomes supersaturated to a certain degree. 
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Fig. 2. — Spatial profiles of (a) temperatures (T g : gas, T^. dust, T ra d: ambient radiation) and 
(b) densities (pd ; dust, p v : silicate vapor) in the vicinity of the shock front for i? p = 100 km, 
n = 10 15 cm -3 , v B — 8 km s -1 , £ = 0.01, and a d = 1 /mi. Panels (c) and (d) are, respectively, 
expansions of panels (a) and (b) in the distance scale. in panel (a) indicates the dust 
temperature at the beginning of the second stage (see text). 
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Fig. 3. — Evaporation fraction rj as a function of dust temperature at the second stage, 
Td2- We take T^2 at the time when the velocity of the dust particles relative to the gas is 
1/10 times thermal velocity of the gas molecules. The symbols indicate numerical results 
and the curves show rj calculated with use of Eq. Q together with Eq. @ for given At, 
cooling timescale of the gas. Plotted in the panels are (a) numerical data for all sets of the 
parameters, (b) those distinguished by the _R p -values given in the panel, (c) those by no, and 
(d) those by £. 
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Fig. 4. — Comparison of T d2 , dust temperature in the second stage, obtained by the nu- 
merical simulations (vertical axis) and those calculated with the use of Eq. (IB2I) given in 
Appendix IB1 (horizontal axis). 
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Fig. 5. — The same as Fig. |3]but the values of in the horizontal axis are replaced by 
those calculated by using Eq. (IB2I) . 
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Fig. 6. — Temporal variations of (a) the radius R of the gas cloud and (b) its expansion 
velocity v r behind a planetesimal bow shock. The curves for 7 = 7/5 correspond to adiabatic 
expansion of a gas composed of H2 molecules. All quantities including time t are normalized 
(see text for details). The solid curves show exact solutions given by Eq. ( 1C2I) for R and 
by Eq. (IClj) for v r , while the dashed curves show approximations in the two limiting cases 
given by Eq. ( !T3|) . 
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Fig. 7. — Cooling rate of the vapor produced by a planetesimal bow shock during its expan- 
sion versus the gas temperature T. The horizontal axis is T normalized by the initial temper- 
ature T and the vertical one is the cooling rate normalized by T /t s0 , where t s0 = Ro/c s0 is 
the sound-crossing time. The solid, dashed, and dotted curves show cooling rates for 7 = 5/3, 
7/5, and 1.1, respectively. The gray region indicates a range of the equilibrium condensation 
temperatures of silicates under the total pressure of protoplanetary disk, T e = 1300 — 1500 K. 
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Fig. 8. — Relation between the cooling parameter A and the evaporation fraction r]. Plotted 
are rj calculated by the numerical simulations and A calculated from Eq. (ID 2ft : (a) for the 
data for all all of the parameter sets, (b) for each R p , (c) for each no, and (d) for each 
£. The gray region (A < 1) indicates the region where condensation through homogeneous 
nucleation does not take place during the vapor cooling. 
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Fig. 9. — Typical radii a m of the particles condensed in the vapor and the degree of the 
supercooling AT/T e versus the cooling parameter A and the dimensionless surface tension V. 
The gray region indicates the ranges of A expected from the model and of possible values of 
T of the condensates. Solid curves indicate and AT calculated based on the homogeneous 
nucleation theory for T = 10 to 60. Dashed lines indicate those for A = 1 to 10 10 . Products 
of the evaporation and condensation experiments are shown by the red (enstatite crystals, 
Yamada 2009) and green regions (forsterite crystals, Kobatake et al., 2008). The typical size 
range of ultra-fine particles in the matrix of Allende meteorite (Toriumi, 1989) is shown by 
the blue region. 
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Fig. 10. — Formation of chondritic materials produced by planetesimal bow shocks in the 
course of planetesimal growth. 



